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Abstract 

The problem of calculating of the mass spectrum of the two-body Bethe- 
Salpeter equation is studied with no reduction to the three-dimensional ("qua- 
sipotential") equation. The method to find the ground state and excited states 
for a channel with any quantum numbers is presented. The problem of the 
confining interaction for the Bethe-Salpeter equation is discussed from the 
point of view of formal properties of the bound state spectrum, but with only 
an inspiration from the QCD. We study the kernel that is non-vanishing at 
large Euclidian intervals, Re — » oo, which is constructed as a special limiting 
case of a sum of the covariant one-boson-exchange kernels. In the coordinate 
space this kernel is just a positive constant and corresponds to the kernel 
oc S(kE) in the momentum space. When the usual attractive interaction is 
added, it is found that this kernel is similar in its effect to the non-relativistic 
potential in coordinate space, V(r), with V(r —> oo) — > V^. The positive 
real constant Voo gives the scale that define the limit of bound state spectrum 
compared to the sum of the constituent masses, M < 2m + V^. At the same 
time the self-energy corrections remove the singularities from the propagators 
of the constituents, i.e. constituents do not propagate as free particles. Com- 
bination of these features of the solutions allows an interpretation of this type 
of interaction as a confining one. The illustrative analytical and numerical 
calculation are presented for a model of massive scalar particles with scalar 
interaction, i.e. the "massive Wick model". 
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1 Introduction 



The Bethe-Salpeter (BS) equation |T], [2j provides a consistent covariant descrip- 
tion for bound state of two fermions (or fermion-antifermion). In particular, a 
phenomenological success has been achieved in the case of the deuteron describing 
ground state properties and various reactions |3|, ||. It has been established that 
BS formalism reproduces the bulk of the non-relativistic approaches at the nucleon 
momenta small comparing to the nucleon mass, |p| < m, but displays a relativistic 
effects at larger momenta |p| > m and in some of the spin- dependent observables. 
However, the relativistic effects are not too large, since the deuteron is essentially a 
non-relativistic system. 

Another topical subject studied within the BS formalism is the bound state of the 
quark- ant iquark systems, mesons (qq). Compared to the deuteron case the gg-bound 
states are governed by a more complex dynamics and, correspondly, physics here is 
richer. Usual way to proceed here is to reduce the BS equation to an approximate 
equation of the quasi-potential form and then to use approaches and methods similar 
to those which are used in the non-relativistic Schrodinger equation. 

A number of the specific questions of principle and technical problems must be 
resolved in order to reach the same state-of-art as it is in the case of the deuteron: 

1. First, we need a rigorous method to obtain the bound state spectrum and the 
corresponding amplitudes for the BS equation. Today we do not find in the 
literature a discussion of this problem without making approximations. 

2. Second, the fundamental problem here is a construction of an effective theory, 
which would provide the BS equation with confinement. Moreover, it is not 
clear how the confinement should manifest itself in this formalism. 

3. Third problem, we would like to mention here, is the "problem of light quarks 
(u and d)" . For the formalism ofthe BS equation we understand this problem 
as a description of the dynamical generation of the masses and creating the 
bound states with a mass, which is much larger then sum of the "bare" masses 
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of the constituents. Apparently, physics of light quarks is essentially different 
of physics of heavy quarks. We will not discuss the special problems of light 
quarks physics in this paper. 

All these problems were discussed starting the pioneering work of Llewellyn- 
Smith ||. There are, of course, methods to study the spectrum based on the 
reduction of the two-dimensional integral equations to an infinite system of one 
dimensional equations by means of the hyperspherical harmonics and, then cutting 
off the system at some lowest functions as an approximation H |[| [7j . This method 
is valid for the tightly bound states, with mass vanishing compared to the sum of 
the constituent masses, in which case the BS equation has an exact 0(4) -symmetry. 
This was sufficient for early studies, where the mechanism of the confinement was 
modeled as the binding of the heavy quarks to the relatively light bound states. 
However, it has become clear that even for heavy quarks the mass spectrum of the 
mesons lies near and even above the sum of the constituent masses. Therefore other 
methods must be employed to at the spectrum with such characteristics. 

Alternative approaches to the study of the spectrum of the relativistic bound 
state is based on the reduction of the BS equation to one of the possible quasi- 
potential equations and then utilizing the Schrodinger-like potentials, V ~ ar, in 
the coordinate space || or its generalization in the momentum space J9], [10, |llj. In 



spite of the obvious advantages of these approaches, such as clear physical interpre- 
tation of interactions and solutions with relationship to the Schrodinger equation 
approach JT2], still there is interest to study the full BS equation, which gives direct 
connection to the covariant field theory. 

As to the underlying effective field theory for the relativistic equations of the 
bound states, now we can only imagine about the structure of such a theory and its 
connection to the fundamental theory, QCD. The approach based on the integral 
equations of the field theory [|13|], the set of Schwinger-Dyson equations for prop- 
agators, the Ward-Takahashi relations (Slavnov- Taylor relations) and the Bethe- 
Salpeter equation is very promissing. The model investigations presented here al- 
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low for a study of the properties of the analytical structure of the quark propaga- 
tor [FTS], [I3J, the mechanisms of the chiral symmetry breaking and the dynamical 
generation of the masses and the connections with the bound state problem |I5| 
and observables of reactions 



In the present paper we (i) consider the problem of the spectrum of the two-body 
BS equation with no reduction to the quasipotential equation or other approxima- 
tions and (ii) attempt to construct a model for the confinement, based on the formal 
properties of the kernels, but with no direct connection with QCD. We present, in 
Section 2, a method to find the ground states and excited states in the channel 
with any quantum numbers. The illustrative numerical calculations, aimed to show 
capabilities of the method are presented in Section 2.3. We consider in Section 3 
the formal properties of the kernel in the form of a sum of the ladder kernels. We 
derive conditions when this kernel displays properties, similar to the non-relativistic 
constant potential in the coordinate space. The implementation of this interaction 
to the single particle propagator leads to the disappearance of the singularities in 
the propagator. Using the developed method to study the bound state spectrum we 
analyze the BS equation with this kernel and self-energy corrections in Section 3.3. 
Combination of properties of the solutions in this case allows for an interpretation 
of the interaction of this type as the confining interaction. 

Nota bene. In this paper we often use the name "quark" for the constituents and 
the "meson" for the bound system. This is just a matter of convenience and our 
concrete calculations are not based on the QCD and do not pretend to a phenomeno- 
logical application. Moreover, the illustrative analytical and numerical calculations 
are presented for a model with massive scalar particles with scalar interactions, 
the "massive Wick model". However, the methods can be applied directly to the 
BS equation for two spinor fields interacting via a phenomenological gg-interaction. 



The spectrum of the spinor-spinor BS equation has been studied in ref. []17j with no 
confining interaction, and the application of the model with the confinement to the 
meson spectrum will be done elsewhere. 
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2 The bound state spectrum of the Bethe-Salpeter 
equation 



2.1 The scalar Bethe-Salpeter equation 

The BS equation for bound state of two equal mass scalar fields, "quarks", interact- 
ing via scalar exchanges has a form: 

^{p ) P) = D{p l ,p 2 y 1 J jULjC(p,p',P)iP(p',P), (1) 

where i[>(p, P) is the BS amplitude for the bound state with momentum P, JC{p, p', P) 
is a kernel of BS equation and D(pi,p 2 )~ l is the two-body propagator. Generally 
speaking, the kernel K, is a sum of the all renormalized irreducible graphs presenting 
the Green's function with two incoming and two outcoming constituent fields and 
the propagator D" 1 is product of the renormalized full single-particle propagators. 
In practice, however, both of them are usually taken in the lowest order in the 
coupling constant. It corresponds to the ladder approximation for the kernel with 
free particle propagators. In the simplest case of an exchange field of mass \x and 
coupling constant g, the ladder approximation for the kernel has an explicit form: 

i g 2 

^ladder {^,P,p',P) = ^ladder 0, P, P') = 1 77^ (2) 

(p — p') 2 — fi 2 

The two-body propagator then reads 

D{pi,p 2 ) = d(pi) • d(p 2 ) = (p\ - m 2 ) • (p\ - m 2 ) , (3) 

where dip)" 1 is free one-particle propagator, m is the mass of the constituents and 
pi t 2 = \P ± p are the constituent momenta. 

Phenomenological applications of the equation (JU) (or more interesting spinor- 
spinor equation) requires a more general form of the kernel than eq. (0). More 
degrees of freedom are provided by introducing a sum of the exchanges with different 
parameters and different Lorents structure, e.g. vector or pseudoscalar exchanges. 
In the presence of the several (effective) exchanges with different quantum numbers 
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and in higher orders kernel may contains both attractive and repulsive terms. Here 
we consider the kernel of the form: 

N ig 2 

where index j enumerates different "exchange" terms, distinguished by the mass 
parameter, p,j, and strength, gj. Factor ej = ±1 defines attractive (e, = +1) or 
repulsive (e/ = — 1) type of the correspondent term. Being interpreted as contribu- 
tion of the lowest order diagram, the terms with negative value of ej require either 
antihermitian term in the interaction hamiltonian (lagrangian) or fields with addi- 
tional quantum numbers. However, we consider the entire sum in the eq. @ as 
a convenient parametrization of the unknown full kernel of the BS equation. The 
convenience of such a form of the kernel is that kernel is explicitly covariant and con- 
tains only "field theoretical" elements, the free particle propagators. Accepting this 
way of action, we assume that parameters of the kernel should be fixed to describe 
the experimentally known spectrum of the system. 

The free propagator D(pi,p2)~ 1 of the form (§) and the kernel in the ladder 
approximation, (fj) or (||), usually provide the basis for an investigation of the BS 
equation. The parameters for a phenomenological adjustment of the theoretical 
framework are the coupling constants, gj, and exchange masses, The parameter 
m in the propagator is referred to as "physical mass" of the constituent and it 
is supposed to include effectively the self-energy corrections. In the case of the 
observable particles, e.g. nucleons, m is the observable mass, otherwise, e.g. quarks, 
it has rather ambiguous model-dependent value. 

In order to study the spectrum of the BS equation we, first, fix the frame of 
calculation as the rest frame of the system, where P = (M, 0) and M is the mass of 



the system, "meson". Then, we perform the well-known Wick rotation fll8l , which 
may later cause difficulties in the calculation of observables of some reactions, but 
does not affect the spectrum of the equation and simplifies the numerical analysis. 
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Under the Wick rotation the BS equation keeps the form (|l|) with the ladder kernel 

f 

(p-p') 2 E + /i 2 ' 

and propagator 



Kladder(jt,P,p') = 7 7^—, o ' (P ~ P')e = (Po ~ Po) 2 + (P ~ P') 2 > ( 5 ) 



D(p, M) = [(p 2 + m 2 + p 2 - M 2 /i) 2 + M 2 p 2 \ . (6) 

We use the previous notation for the "rotated" functions ip, JC and D, and for the 
"Euclidean" momenta, p and p' . 

Next, we perform a partial wave decomposition of the equation P"D| |: 

1 2tt 00 L 

(v ^2—772 = | D . i d ,i E E QL(f3)YLx(9p,V P )YL(9> p ,v' p ), (7) 

\P-P )E+ V I P I ' I P I L=0 X=-L 

1 oo £ 

^(p,M) = n £ E ^L(Po,ipi,M)F iA (^,^), (8) 

I P I L=0\=-L 

where we have already taken into account that ipL are independent on the projection 
of angular momentum, A; the dimensionless parameter j3 defined by the expression 

fi 2 + p 2 + p' 2 + (go - p' f 

I o I I I /I ' 

2 I P I • I P I 

and are the Legendre functions of the second kind. For I = 0, 1 they are: 

= i.n g±l) , ftM = |ln(j±l)-l. (9) 

Substituting ([7]) and (§) into eq. (|T|) and performing the angular integration, we 
arrive at a set of the independent equations for the partial amplitudes ipi, corre- 
sponding to the states with the angular momentum L: 

MPo, |p|, M) = D(p, m)- 1 J ^0±JC L (p o , |pU, \p'\)Mp'o, Ip'I, m),(io) 

IC L (Po,\p\,PoAp'\)=9 2 Ql((3). (11) 

The amplitude of pure quantum state with defined angular momentum L and its 
projection A is then given by 

V(p, M, L, A) = Vx(Po, |p'|, M)F LA (^, P ). (12) 
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The two-dimensional integral equations ( |T0| ) are an exact projection of the initial 
equation (p]) with @ and (D). In the case of the generalized kernel @, there will 
be a sum over different terms on the right hand side of eq. (jTD|): 

K%(po, \p\,p' , |p'|) = J2e J g*Q L W J ). (13) 

i=i 

To obtain the spectrum of the BS equation we need to solve the eigenvalue 
problem for the bound state mass M at fixed set of parameters of the kernel. =It is 
to be noted that the BS equation is not linear in the mass M. On the other hand it 
is linear in the coupling constants squared, g 2 . Therefore we consider an equation 
of the form: 

MPo, |p|,M) = A • D(p, M)- 1 J ^M/c§(pd, IpU, Ip'I)^(Po, |P'I,M), (14) 
or symbolically 

i/j = X-JCip, (15) 

where we skip indicies L for shortness. 

Solving the linear eigenvalue problem for A at fixed parameters of the kernel and 
various values of the bound state mass, we can map A(M). Then by inverting the 
mapping as M(A), we will find solutions of the eigenvalue problem for M with the 
kernel A-/C, where factor A is trivially absorbed by redefining the coupling constants. 

Traditional methods to find solutions to the linear eigenvalue problem corre- 
sponding to the integral equation ( |l"5"D are based on the ideas of the classical Fred- 
holm theory. The basic idea here is to substitute integration by a summation and 
then deal with a sufficiently large system of linear algebraic equations. The prob- 
lem is reduced to finding the eigenvalues of the matrix corresponding to the integral 
operator /C. However, in the case of the covariant BS equation such a program of ac- 
tion meets extra difficulties in view of the double integration, which leads necessarily 
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to very large matrices, say ~ 10 x 10 4 or even larger in the case of spinor-spinor 
equation (with additional factor of 64 if the parity is conserved by interactions and 
factor 256 if it does not). 

We use an alternate method based on the iteration of the BS equation in the 



form (|i0l). This method in the essence is similar to the Malfliet-Tjon method |20 



employed to solve the integral Schrodinger equation in the nuclear physics (see also 



discussion in ref. pT 



2.2 The iteration method 



The iteration of any trial function, with eq. (]TD|), is understood as obtaining 
other function $^ n - ) using the algorithm (we use the symbolic notations as in eq. 



3)): 



A-/C 



$(0)_ 



(16) 
(17) 



In order to organize the iteration (|i~6|)- (|17D on the computer, we need an "in- 
tegrator", corresponding to the operator JC. This means that we have to perform 
the computer calculation of the double integral on the r.h.s. of the equation with 
a defined kernel and any trial function, (we assume good enough behavior 

of the function, such as absence of singularities and vanishing at large arguments, 
|p| — > oo and po — > ±oo). In our particular calculations this integrator is organized 
as a two dimensional Gauss integration with suitable mapping of variables. Next, 
we would like to know what happens with equation after sufficiently large number, 
n, of iterations. 

Let us assume that solutions, ipa, of the equation ( |15"D corresponding to the 
eigenvalues X a belong to the complete system of functionsf22|. Therefore, the trial 
function can be expanded as 

oo 

$(°> = $> a Vv (18) 

Q = 
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Thus, the result of iteration (pf)-( |l7|) is: 



a=0 



A, 



"00 



From last equation it is obvious that at sufficiently large n 
a > are small compared to the ground state term, a = 0: 



lim $ (n) = $ (A ° = C ■ ip + O 

n->N 



Ao 
Ai 



(19) 

iV all terms with 
(20) 



Therefore N to be chosen to make the last term on the r.h.s. of to be negligibly 
small. Then, comparing $W and c^^ 1 ), we find the ground state eigenvalue, Aq: 



Ao 



$(JV+1) 



(21) 



This recipe to find the ground state eigenfunction and eigenvalue works nicely 
numerically and, of course, does not depend on the choice of the initial trial function. 

Formula flT8|) also provides us with the possibility to find higher levels on A. 
Indeed, taking the combination of iterating functions: 



A_ 
A7 



. $» = x; ( i - 



a=l 



A a 

a7 



A, 



A,, 



(i+l) 



(22) 



i.e. new trial function which does not contain admixture of the ground state! It 
is easy to see that iterations of this function give the eigenfunction ipi of the first 
excited state, Ai, similarly to the procedure for ground state. The same procedure, 
in principle, can be organized for any desired level on A. Thus, the problem is solved. 
The only limitation is, of course, accuracy of the numeric calculations. Calculation of 
the high A require a precise calculation of all levels below, which leads to substantial 
computer time. No special numerical problems were find in calculating the three 
lowest levels for eq. ([15]). The numerical results are presented in the next section. 



11 



2.3 Numerical results for iteration method 

To study numerically the capability of the method to solve the integral BS equation 
we consider eq. (|T0|) with the model kernel of the form (|T3|). The parameters of the 
model are presented in Table I. We refer to the constituent fields as "quarks" and 
their mass, m, is chosen to be similar to the mass of the c-quark. The parameters 
of the kernel are chosen to provide the typical density of the levels of the cc-bound 
state, the charmonium, not too far from the limit of the spectrum Mi im = 2m. 



coupling constants 


e i 


mass 


g]/(4n), GeV 2 




fij, GeV 


37800.0 


+1 


0.10 


37800.0 


-1 


0.11 


45.0 


-1 


0.95 


45.0 


+ 1 


1.425 


m = 1.5 GeV, 



Table. I. The parameters of model ("scalar charmonium"). 

Results of a calculation of the spectrum, A(M), are presented in Fig. 1 for the 
three lowest levels in the channels with L = (S-states) and L = 1 (P-states). The 
physical spectrum M(A) can be obtained crossing the plot by the line A = Ap^. 
For instance, as is shown in Fig. 1 at \ p h ys = 0.13 the masses of the lowest states 
are (with accuracy ~ 0.5%): 

M(1S) = 1.265 GeV; M(2S) = 1.939 GeV; M(3S) = 2.251 GeV; (23) 
M(LP) = 1.751 GeV; M(2P) = 2.146 GeV; M(3P) = 2.385 GeV. (24) 

For illustartion, we pesent the amplitudes ipo, corresponding the spectrum (p3|) 
in Fig. 2. These amplitudes are shown as a functions of the spatial momentum |p| 
and at po — 0. We see that the type of radial excitations is similar to the one for 
the Schrodinger equation. However, in general case of the BS equation we deal with 
the hyperradial excitations of two-dimensional surface, iPl(Po, |p|)- 



12 



3 Confinement for the Bethe-Salpeter equation 



3.1 General discussion and non-relativistic confining poten- 
tials 

The idea of confinement has different realization within different theoretical ap- 
proaches. The simplest intuitive picture is given by the non-relativistic bound state 
formalism based on the Schrodinger equation with a QCD inspared phenomenologi- 
cal potential. A system of two particles interacting in a non-confining potential, van- 
ishing as r — > oo, has a spectrum of bound states with an upper limit, Mn m = 2m q , 
and a continuum above this limit. The confinement is conventionally associated with 
the infinitely rising linear part of the full gg-potential, V\ = ar, a > 0, which pro- 
vides with the mass spectrum extending infinitely beyond Mi im . It is clear that this 
mechanism can not be directly adopted by the covariant field theoretical approaches, 
such as the BS formalism. 

More relevant approach is based on the simultaneous analysis of the BS, Schwin- 
ger-Dyson (SD) equations and Ward-Takahashi at the lowest order and with the 
model gluon propagator [|15|, In particular, the role of the analytical structure 



(structure of singularities) of the quark propagators is discussed here. It is found 
that, with certain choices of the model gluon propagator, the quark propagator is an 
entire function (function with no singularities) at physical momentum of the quark. 
This important property of the quark propagator is considered as an indication of 



confinement O, 04 16 



These two examples present two essentially different pictures of what is referred 
to as confinement. In the first case, the confinement is the two-body effect, i.e. 
quarks can not escape from each other because of the interaction between them. 
In the second example, confinement is attributed to the property of a single quark, 
which can not propagate clS db free particle. 

The general formalism of the field theory suggests that the two-body and the 
single particle phenomena are in a generic relationship. So do the approaches based 



13 



on the BS and SD equations. Our approach to a modeling of the confinement is in 
some sense inverse to that of refs. fllBI, HM. We, first, construct the covariant kernel 



of the BS equation which, we expect, would provide the confinement and only then 
study the modifications of the quark propagators involved in the BS equation. 

We start from a few unsophisticated observations prompted by the non-relativistic 
picture: 

1. In the non-relativistic limit the covariant theory can be reduced to the for- 
malism with the Schrodinger equation, where we can expect the picture of the 
confinement as interaction with non- vanishing potential at r —>■ oo to be valid. 
The main distinguishing feature of the spectrum here is the existence of the 
bound states above the two quark mass limit, M Hm . 

2. It is not necessary to have an infinitely rising potential, if we intend to discuss 
only a few lowest levels in the spectrum. More manageable potential, V — ► 
Voo > at r — > oo, could be sufficient, if is large enough. In this sense, we 
also refer to the constant potential as a confining one. 

3. The potential in the momentum space V(k) can be obtained as a non-relativistic 
limit of the kernel /C, similar to 

nk) = (25) 

if the kernel is of the form @. Important point here is that the non-relativistic 
form of the potential (EH) is of a field-theoretic origin. 



Basing on the last observation we expect that if we define a way to construct a 
confining potential from the non-relativistic field theory, then we can apply similar 
methods to obtain the relativistic confining kernel. Very often the following recipe 
is used. The non-relativistic Yukawa potential in the coordinate space V(r), the 
Fourier transform of eq. ( f25|) , is 

a 2 e~^ r 

V{v) = -f- — . (26) 
4tt r 
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The linear potential can be derived as 

d 2 



V, 



lim 



dfj,'- 



-V(v) 



lim — re Mr 

fi^O 47T 



47T 



-r. 



(27) 



The the relativistic generalization is made by a Fourier transform to the momentum 



space and replacing the non-relativistic Yukawa potential, fl25), by it relativistic 
analog: 



K oc lim 



<9 2 



<9/i 2 A; 2 — /x 2 



(2* 



Taking the limit in (p8|), one should exercise great deal of care, since this leads to 
the appearance of generalized functions in the kernel || [TT|j . 

This recipe give us a guideline, however it is not completely satisfactory, since 
(i) the kernel (|28|) (or the potential (^7|)) is not of the field-theoretic form and (ii) 
it is not clear does the direct use of operation (^J) lead to the rising or, at least, 
non- vanishing interaction in the four-dimensional space. 

Intending to stay with our parametrization of the kernel as a superposition of 
the ladder terms, similar to (§), we have to find an appropriate presentation of the 
operation (^). Let us start with a superposition of the non-relativistic potentials: 



V{v) 



E y ex p[ 



-fxajr\ 



(29) 



where fi provides the mass scale and aj are dimensionless parameters. Then ex- 
panding the exponents we get 



V(v) 



1 2 a 2 r 

- - w + v ~y 



(30) 



From eq. (|30| ) we see that, taking the limit fx — > and correspondly adjusting 
the parameters Cj and aj, desired non- vanishing behavior of the potential can be 
provided. For instance, for constant potential, V c , we have 



(31) 



15 



where we need only two terms to satisfy the conditions, i.e. j max = 2. For linear 
potential,!^., 

c i = /x- 2 Q ; EQ = °; E^4 = °; E a ^ = ^^ ' ( 32 ) 
j j j 

where j max = 3. The limit fi — > corresponds to the physical picture of the superpo- 
sition of very light mass exchanges with slightly different masses and large coupling 
constants. Please, note that power of the non-vanishing term in the limit ft — > is 
solely controlled by the power of the fi in the denominator of the coefficients Cj. 

Using a superposition of the Yukawa potentials in the momentum space, we find 
(see Appendix A): 

V c (k) = A c 6®(k), (33) 
Vm = (I_^. V ,, (k) _|. (34 , 

Note that, for the linear potential, the Schrodinger equation becomes an integro- 
differential equation in momentum space. 



3.2 Confining kernel for the Bethe-Salpeter equation 

We look for a confining kernel, /C con , of the BS equation in the form of the 
superposition of ladder kernels in the momentum space and after the Wick rotation: 

C ■ 

where k% = k\ + k 2 . 

In spite of the obvious similarity to the non-relativistic case, the expression for 
the relativistic kernel, fl3"5p, has essentially different properties in view of the larger 
dimension of the space. Indeed, the ladder kernel in the coordinate space is [|J: 

\r (t> \ 2 Ki(/xi? s ) 

>^ladder{-tlE, A 4 J — 9 f^ 5 , [oV) 

Ke 
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where Re = (r 2 + t 2 ) 1 / 2 . The asymptote at small Re can be obtained as: 

K-ladder{RE ~ » 0, /i) ~ g 2 I 

[ K E 

4± ^TT^ h (^) - 5 <« m + '> + « m + 2 »1 } ' < 37 » 

where ^(m) is the Euler's psi function. The expansion ( p7[ ) is different from the 
behavior of the non-relativistic potentials. However, a procedure similar to that 
of the non-relativistic case can be applied to cancel the lowest order terms in the 
expansion (|37D and generate a non- vanishing kernel in the limit [i — > 0, since at large 
Re there is the exponential suppression similar to the one in the Yukawa potential: 

/Cladder{RE ~> OO, /i) ~ g /i 1 ^2 ' ( 38 ) 

Re 

We can see that the direct use of operation fl28|) does not lead to the desired non- 
vanishing at Re — ► oo behavior of the kernel. 

We study here the kernel with the lowest power in which as we expect 
controls the asymptotic behavior at large Re- Analysis in the momentum space, 
similar to that of the non-relativistic case, gives the following conditions for the 
coefficients Cj (see Appendix B): 

c j = ^Cf, Cj = °; E a Pj = °; E «| |na ^ = A ^ °> ( 39 ) 

J 3 3 

which can be satisfied explicitly for j max = 3. Note, that the second condition 
cancels the most singular terms, ~ R E 2 -, in the coordinate space, fl37|), similar to the 
non-relativistic potentials. This choice of conditions leads to a kernel of the form: 

JC con {k) = -{2ixfU^\k E ), (40) 

where for simplicity we introduce new effective coupling constant U, which has di- 
mension of mass, and sign is chosen to provide us with a confining- like kernel, similar 
to the positive constant potential in the non-relativistic case. Fourier transform of 
the kernel, (|40D, is a constant, ~ f/ 4 , in four dimension. This means that it does 
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not behave like a non-relativistic constant potential, for which we expect behavior 
like ~ S(to) ■ constant. Therefore, the kernel fl4"0| ) does not exactly correspond to 
the non-relativistic constant potential and the effective constant, U 4 , is not related 
to constant in such a potential. Note, since the constant f/ 4 does not have the di- 
rect physical meaning, the choice of the factor (27r) 4 is arbitrary and it is made for 
further simplification of formulae. 



The form of the kernel, (40), in accordance with our main idea, is considered as a 
special limiting case of the sum of the ladder kernels (sum of one-boson exchanges), 
which provides explicit covariance of the kernel and connection with the usual field- 
theoretic constructions. (Note that this is a valid form in the Euclidian space, 
whereas the transition to the Minkowsky space is not defined.) By itself the 5-form 
of the kernel is not something very unusual in studying of the bound states of the 
quarks. For instance, such a form is considered as "regularized" form of the highly 
singular, ~ k~ A , behavior of the gluon propagator [O and is a basis of the models 



for studying of the SD equation for the quark propagator |L5| In particular, this 
form of the gluon propagator leads to the quark propagator without singularities 
along physical momentum. In the lowest order, such a gluon propagator gives the 
kernel of the BS equation JTH . 

Let us study the effect of the kernel, (|40|) , in the BS equation, ([I]), under the 
Wick rotation. The form of K, con allows an integration in the equation explicitly: 

^• p) = -i^) iip ' n (41) 

From eq.pT|) we find that the kernel, JC con , does not allow for the bound state 
solutions of the BS equation. In this sense we can expect that this kernel, in effect, 
is similar to the constant non-relativistic potential, which along does not allows for 
the bound states. If this is the case, we expect that adding this kernel to the regular 
attractive kernel, i.e. like the one presented in Section 2.3, Table 1, we will get a 
shift of the spectrum by a constant. 

The BS equation with the combined kernel, Kq + )C con , can be transformed to 
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the usual form, but with a modified two-body propagator: 

"'("• P) = D(J 2) + u* JwF ,CG(p ' p '' PHip '' n (42) 

New propagator in eq. (E5J) has different analytical properties, compared to the 
initial free propagator under the Wick rotation, (^]), and main difference is that the 
new propagator does not have singularities at M > 2m, which were a signal of the 
limit of the spectrum at this point. This is an indication that the physical spectrum 
exists beyond this point. However, without numerical analysis we are not able to 
discuss the properties of the solutions at M > 2m. 

Before to go to the numerics, let us discuss the possible effects of the self energy 
corrections in the presence of the interaction of the form (^0|). Dealing with the 
kernel (HOT) as a sum of the lowest order (ladder) kernels, we calculate the one-loop 
self-energy corrections to the single quark propagator, d{p)~ 1 . Integration over the 
loop is performed at imaginary po component of the four-momentum of the quark 
and result can be analytically continued to any values of Pq. For physical momentum 
we get: 

d(p) = p 2 - m 2 + -. (43) 

The propagator (|43|) does not have singularities for physical values of the momen- 
tum, p. This property of the interaction of the form ([It]) has already been established 
within model investigations of the behavior of the quark propagator [[14], [13). Now 
we see that this single quark effect is in a generic relation with the two-body confin- 
ing interaction in the framework of the BS equation. For our calculations it is also 
important that the singularities of the modified two-body propagator in ( |42] ) with 
D{piiP-i) = d{pi) ■ d{p 2 ) defined by eq. ([43]) still allow to perform the Wick rotation. 



3.3 Numerical investigation of the Bethe-Salpeter equation 
with the confining kernel 

First we compare the meson spectra obtained with eq. (f42|) with U 4 = and U A ^ 0. 
For convenience we take the same kernel, JCq, as that in the section 2.3. The 
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main effect, we expect here, is the shift of the spectrum beyond the limiting point 
Mu m = 2m. In order to compare the spectra, we calculate A(M) for two lowest 
states in the channels with L = 0, 1. This is enough to make conclusions about (i) 
the position of the bound states, (ii) limiting point of the spectra and (iii) separation 
between ground and excited states. Results of these calculations are presented in 
the Fig. 3, groups of curves A and B. The constant U for the case B is chosen to be 
of the typical energy scale of the equation, U = m. 

We find, indeed, that the masen spectrum of the equation with U ^ is extended 
beyond the point Mu m . However, it displays unusual behavior beyond the point, 
M^ m (this point is shown on Fig. 3 by the arrow) |2"3fl . Non-monotonical behavior of 
the curves A(M) indicates some difficulties in the interpretation of the correspond- 
ing solutions. The obvious difficulty is the existence of two solutions with the same 



coupling constant A and different masses, M, (see e.g. discussions in refs. [19 



Therefore, we have to find a way to isolate the only one physical solutions. The so- 
lution of this problem is quite simple. Calculating the normalization of the solutions 
of the equation along M we find that the ground state solutions beyond Mj? m have a 
negative norm, i.e. are the abnormal non-physical solutions. f24|. This observation 
gives us the selection rule to eliminate the extra solutions. 

Another problem with the solutions corresponding U A ^ is that spectrum is 
not only shifted to larger masses, but also the separation between levels is drastically 
increased (see Table II). This effect make it a problem for a phenomenological use 
of the kernel of the form (fTO|). Indeed, if we intend to consider states above Mu m 
we have to take U A large enough to provide us with a new limit of the spectrum, 
however this can be in conflict with the desired separation between levels. One 
may try to adjust the parameters of the remaining part of the full kernel, JCq, so 
as to have a reasonable density of levels. However our analysis showed that in the 
presence of the kernel, JC con , separation between levels depends only weakly on the 
kernel, K,q- 
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A M 


A 

^phys = 0- 13 


B 

\hys = 0- 25 


C 

^ hys = 0-44 


D 

^hys = 1M 


IP - IS 


0.49 GeV 


0.77 GeV 


0.48 GeV 


0.48 GeV 


2S - IP 


0.19 GeV 


0.47 GeV 


0.16 GeV 


0.17 GeV 


2P - 2S 


0.21 GeV 


0.42 GeV 


0.22 GeV 


0.26 GeV 



Table. II. The levels splitting in the "scalar charmonium" for different kernels . 

It should be remembered that in the BS equation the addition of interaction of 
any kind to the full kernel cannot be related linearly to the shift in the mass of the 
system. The same interaction also "shifts" the masses of the constituents through 
the self-energy corrections to the single-particle propagators. A self-consitent ap- 
proach has to be adopted to include both the self-energy corrections and changes in 
the two-body interactions. 

We take into account the self-energy correction, (f43l) , to the quark propagator. 
This leads to a corresponding modification of the two-body propagator D(p, 
(|3|). The resulting spectrum is presented in Fig. 3, group of curves C. It is clear that 
for the interval of masses ~ 2 4- 3.5 GeV the density of levels "returns to normal", 
of the same value as in the case A (see Table II). The value of constant \p hys = 0.44 
is chosen as an example giving the spectrum density close to the the one of the case 
A. At smaller Xp hys the separation of levels is even smaller. 

The examples of calculations, B and C, show that the kernel containing the part 
)C con is indeed similar in its effect to the non-relativistic potential in coordinate space, 
V(r), with V(r — > oo) — ► V^. where the positive real constant defines the shift 
of the bound state spectrum compared to the case with no K, con , the case A. From 
Fig. 3 another similarity with the non-relativistic constant potential is obvious, this 
kernel gives only a limited number of states in the spectrum, which can be adjusted 
by varying the \ p h ys - For instance, for the case C with Xp hys = 0.2 there is only one 
bound state, IS, whereas with Xp hys = 0.45 there are four states in the S and P 
channels (there can be other undetected states in channels with higher L — 2, . . .). 

On the other hand, taking account of the self-energy corrections, corresponding 
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to this type of interaction, leads to the disappearance of the poles in the single quark 
propagator. Therefore quarks cannot propagate other than being bound in a bound 
state. Apparently, this fact is not related to the infinitely rising interaction between 
them, but rather to the modification of the single quark properties. 

To be sure that the picture is valid in a wide interval of the constant £/ 4 , also 
the spectrum is calculated for the case U 4 = 4m 4 and self-energy corrections taken 
into account. The result is shown in the Fig. 3, group of curves D. That the picture 
is found to be similar to the one of case C, but the spectrum is shifted up to even 
larger Mj? m ~ 4 GeV. The density of levels in the interval of masses ~ 2.5^-4 GeV is 
the same as those of cases A and C (see Table II). The value of constant Xp hys = 0.44 
is chosen again to give the spectrum density close to the the one of the case A. 

Note that in our model the constant U can not be taken arbitrary large, since 
all the calculations are performed in the lowest order. The natural criteria on the 
maximum value of U is that the corrections of the lowest order, ~ ?7 4 , must not be 
too large. For instance, a shift of the mass spectrum with U ^ should not be too 
big compare to the typical masses for the case with U — 0. 

4 Conclusions 

We have presented a method to find the ground state and excited states of the 
two-body Bethe-Salpeter equation for a channel with any quantum numbers. This 
method allows us to solve the bound state problem without reduction of equation 
to quasipotential form or any other approximations. 

Based on a qualitative analogy with the construction of a non-relativistic poten- 
tial with non- vanishing asymptote at large distances, r — > oo, we have proposed a 
recipe to obtain the confining kernel for the Bethe-Salpeter equation, parametrized 
in the form of a special limiting case of a superposition of ladder kernels. We find 
that in the simplest case such a kernel is proportional to in the Euclidean 

momentum space, which corresponds to a constant kernel in the coordinate space. 

We have studied the effect of this kernel on the spectrum of the Bethe-Salpeter 
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equation, when the usual attractive interaction is added. It is found that this kernel 
is similar in its effect to the non-relativistic potential in coordinate space, V(r), with 
V(r — ► oo) — > Voq. The positive real constant gives the scale that defines the limit 
of the bound state spectrum compared to the sum of the constituent masses, M < 
Zm + Voo. At the same time the self-energy corrections remove the singularities from 
the propagators of the constituents, i.e. constituents do not propagate. Combination 
of these features of the solutions allows an interpretation of this type of interaction 
as a confining interaction. 

The illustrative analytical and numerical calculation are presented for a model 
of massive scalar particles with scalar interaction and do not pretend to a phe- 
nomenological application. However, the developed formalism can be straightfor- 
warly adopted for the Bethe-Salpeter equation for the bound state of two spinor 
fields and, therefore, can be used for the realistic studies of the properties of the 
quark-antiquark systems, mesons. 
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Appendix A. The non-relativistic confining po- 
tentials in the momentum space. 

In order to establish the form of the potentials defining by the eq. (p9|), ( |3T| ) and 
( |32|) , let us consider auxiliary integral, I a \ 





I, 



a 








(A.l) 



(A.2) 
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where /(k) is any function for which we assume a "good" behavior (/—►() when 
k — > oo, no singularities, existence of derivatives, etc.), a — c,l depending on which 
of conditions the fl3lf) or (|32|) is imposed on the potential and V a (k) defined as a 
fourier transform of eq. (E^): 



A; 2 + a 2 /i 2 ' 

The limit \x — > is assumed and we take it later. 

Using the common condition Cj = 0, we rewrite eq. ( A.2 ) as 



(A.3) 



dk 



(2tt) 3 ^ ~ J ~> J k 2 + a 2 fi 2 



/(*,«)■ 



(A.4) 



oo 




y cf& arctg 










f(k,n), 



(A.5) 



Integrating by parts the last integral in (|A.4 ), we get: 
where /' = df /dk. 

Let us now consider integration over k only. These integrals on the r.h.s. of eq. 
(|A.5Q can be split in two parts: 

otjfi oo , 

k 



I a oc fj, C j a j If + J \ dk arct S 



f(k,n). 



(A.6) 



3 I ctjii 

We estimate the first integral by the mean value theorem: 



otjfi, 



dk. 



A* 2 (^yli:^^^), 



(A.7) 



J '3 

where < £ < 1. 

The second integral is estimated using the expansion arctgx = tt/2 — 1/x + 
l/(3x 3 ) — . . ., which is valid at x > 1. It can be shown that the first term in this 
expansion gives the leading contribution to the full integral in the limit // — > 0: 



oo 



(At 



dsn 
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Finally, taking the limit /x — > and accounting for the conditions on the coeffi- 
cients Cj, we find 

Ic = AJ(0), (A.9) 

* = Q-v)^ //(0) - (A ' 10) 

These equations give us the potentials in the form (|33|) and fl3~4|). 

Appendix B. The relativist ic confining kernel in 
the momentum space. 

In order to establish the form of the kernel defined by eq. (^) at the limit /x — > 
and lowest (but non-zero) degree of /x -1 in the Cj, let us consider the auxiliary 
integral, J: 

d A k 



( .,_ )4 £co»(A0/(A0 (B.l) 
Jj^J dkk 3 JC con (k)f(k, fi (4) ), (B.2) 



where k is four-momentum in Euclidian space, k = (k 2 -!-/^) 1 / 2 , f^ 4 ) is the hyperangle 
defining the orientation of the vector k in the four dimensional space; f(k) is an 
arbitrary function for which we assume "good" behavior (/ — > when k — » oo, no 
singularities, existence of derivatives, etc.). The limit /x — > is assumed and we take 
it later. We are omitting all nonessential factors, such as 2n, etc., in the following 
calculations. 

Adding and subtracting expression Cj/k 2 to each term in (35|), we rewrite eq. 
(ggD as 

/ oc -/x 2 / ^ (4) E Iwr^rJ^ fi(4) ) 

i o J 

+jd^j:cj^f(k,^). 

3 

Th last term is cancelled by the condition 

E C,- = 0- (B.4) 
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Let us now consider integration over k only. This integral in r.h.s. of eq. (|B.3| ) 
can be evaluated as: 




(B.5) 



since the function remaining under integration has sharp maximum at k — otj/J,. 
We also introduce the cut-off parameter A to regularize formally the logarithmically 
divergent integrals. At a later stage of calculation we take the limit A — > oo, but 
for the moment it is enough to assume A ^> fi. Performing the integration we get 

j 

Applying the condition 

Y, C 3 a) = 0, Yl Cjojlnay ^ 0, (B.7) 

3 3 

and taking the limit \i — > 0, we get 

I oc fi 2 f(0) Y <" ; n jlnn r (B.8) 

3 

which proves the form of the kernel as eq. (f40|). 
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Figure captions 



Figure 1. Mapping of the bound state spectrum of the Bethe-Salpeter equation. 
The physical states correspond to the cutting of the trajectories A(M) by the line 
A = Xphy S - The three lowest trajectories for L = (solid curves) and L — 1 (dashed 
curves) are presented for the parameters of the kernel from Table I. An example of 
physical cut is shown for \ p h ys = 0.13 and corresponding masses are presented in 



dD-(H. 



Figure 2. The amplitudes corresponding to the three lowest S— states (L = 0) 
for the Bethe-Salpeter equation with interaction from the Table I. These amplitudes 
corresponding to the states with masses eq. (|23|) and X p hy S = 0.13. Curves: ground 
state (IS) - solid; first excited state (2S) - dashed; second state (2S) - dotted. 

Figure 3. Mapping of the bound state spectrum of the Bethe-Salpeter equation 
with different kernels. The group of curves, A, is the two lowest state from the Fig. 1 
(see caption). Solid curves present the S-states and the dashed curves presents the 
P-states. The groups B, C and D are trajectories for the same kernel as the group 
A plus the confining part, JC con : B and C with U A = m 4 ; D with U A = 4m 4 . The 
calculations of groups C and D also take into account the self-energy corrections. 
Arrows along the M— axis show the limits of the corresponding physical spectra. 
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